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Abstract 

Confirming a conjecture of Lyng-Raoofi-Texier-Zumbrun, we show that stability of 
strong detonation waves in the ZND, or small- viscosity, limit is equivalent to stability 
of the limiting ZND detonation together with stability of the viscous profile associ- 
ated with the component Neumann shock. More, on bounded frequencies the nonstable 
eigenvalues of the viscous detonation wave converge to those of the limiting ZND detona- 
tion, while on frequencies of order one over viscosity, they converge to one over viscosity 
times those of the associated viscous Neumann shock. This yields immediately a number 
of examples of instability and Hopf bifurcation of reacting Navier-Stokes detonations 
through the extensive numerical studies of ZND stability in the detonation literature. 
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1 Introduction 

In one-dimensional, Lagrangian coordinates, the reactive Navier-Stokes (rNS) equations 
modeling reacting flow for a one-step reaction may be written in abstract form as 

u t + f(u) x = e{B{u)u x ) x + kqip(u)z, 
z t = e(C(u,z)z x ) x - k(p(u)z, 



where u,f,q£ W 1 , B G R nxn , z, k,C,^e R 1 , and k, e > ED |LyZT| [LyZ2"l |LRTZ| [TZ4] . 



Here, u comprises the gas-dynamical variables of specific volume, particle velocity, and 
total energy; z measures mass fraction of unburned reactant or, more generally, "progress" 
of a single reaction involving multiple reactants |FDl |LyZl] ; (p(u) is an "ignition function", 
monotone increasing in temperature, and usually assumed for fixed density to be zero 
below a certain ignition temperature and positive above; q comprises quantities produced 
in reaction, in particular heat released; and k corresponds to reaction rate. Coefficients B 
and C model transport effects of, respectively, viscosity and heat conduction, and species 
diffusion, and e measures relative size of transport vs. reaction coefficients, typically quite 
small. 

A right-going viscous strong detonation wave is a smoth traveling-wave solution 
(1.2) (u, z)(x, t) = (u, z)(x — st), lim (n, z)(x) = (u±, z±) 



2 



of solutions of (jl.ip with speed s > connecting a burned state on the left to an unburned 
state on the right, 



(1.3) z_ = 0, z + = 1, 
with necessarily 

(1.4) <f>(u-) > 0, 0(u+) = = 0, 
and satisfying the extreme Lax characteristic conditions 

(1.5) < • • • < a~_ x < s < a~, a\ < ■ ■ ■ < a+ < s, 

where denote the eigenvalues of df(u±), ordered by increasing real part. Left-going 
viscous detonation waves satisfy symmetric conditions obtained by reflection, x — > —x. 

Multi-step reactions may be modeled by the same equations with vectorial reaction 
variable z £ M m , and coefficients q, C, ip, k modified accordingly. Likewise, the functions 
/, eft, B may be modified to depend, more realistically, also on z, reflecting the different 
chemical makeup of the gas after reaction, with no essential change at a mathematical level. 
For further discussion, see, e.g., (CFj EE EE1 [ZH [LyZTj [LyZ2l ILRTZ1 ITMlHuZ2] . 



A standard simplification in detonation theory is to neglect the small constant e and 
consider instead the formal e = limit, or Zeldovich-von Neumann-Doering (ZND) model 

^ u t + f(u) x = kqip(u)z, 

zt = —k(p(u)z. 

Indeed, there is by now a tremendous body of literature on this model; see for example 
[CFl iFDl iFWl lErTl IEr2l iLSl IBME1 IAT1 lATTl IKS] and references therein. The corresponding 
object to a viscous detonation wave for the (ZND) model is a right-going strong ZND detona- 
tion u° of form (jl.2p - (jl.4p satisfying (II. 6ft . smooth except at a single shock discontinuity at 
(without loss of generality) x = 0, known as a "Neumann shock" \CF\ [Ml EH1 |LyZlj |LyZ2| , 
where u jumps from to u+ as x crosses zero from left to right, with z = 1. We have 
the intuitive picture [CF] of a shock, or "reaction spike" , compressing a quiescent mixture 
and heating it to ignition point, followed by a slow "reaction tail" in which the reaction 
proceeds until all reactant is burned, while, meanwhile, u varies from to u_. 

A ZND detonation profile is determined implicitly \CF\ IHuZ2| by the property, obtained 
by integrating the traveling-wave ODE 

^ -su + f(u)' = kqip(u)z, 

—sz' = —k(f(u)z 

and adding q times the second equation to the first, that — s(u + qz) + f(u) = constant, 
which, together with z = 1 for x > and z{— oo) = implies that 

(1.8) - su + + /(«+) = -su* + /(it*) = -s(u- - q) + /(«-), 
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giving a unique -u_ and profile u , x < 0, for each Neumann shock (u*,u + ) of speed s, so 
long as df(u) — si remains invertible for all < z < 1 along the curve determined by 

(1.9) - s(u + qz) + f(u) = -s(u + ) + /(«+). 

For, solving (jl.9p for u = u(z) by the Implicit Function Theorem then yields the profile 
on x < by solution of the second equation z' = (—k/s)<j)(u(z))z, a scalar equation with 
nonvanishing righthand side, so long as u remains in the region for which 4>(u) > 0. 

A natural question is the relation between the formally limiting (ZND) equations and 
the behavior of the full (rNS) equations as e — > 0. At the level of existence of detonation 
profiles, this was investigated by Majda [M] for a simplified model with u, z £ M 1 and C = 
using direct, planar phase portrait analysis, and extended to the physical (rNS) model by 
Gardner [G] using Conley index techniques and Gasser-Szmolyan [GS] by geometric singular 
perturbation theory. More recently, Williams j W] has revisited the existence problem using 
more quantitative singular perturbation methods, generating detailed matched asymptotic 
expansions to all orders. In each case, with varying levels of detail, the result is that for 
each strong detonation profile (u°,z°) of the (ZND) model with physical choice of /, there 
exists a family (u £ , z e ) of strong detonation profiles converging away from x = as e — > 
to (u°, z°), and near x = to a viscous shock profile for the associated Neumann shock, in 
microscopic variables x = x/e. 

In the present paper, using singular perturbation/asymptotic Evans function techniques 
developed in [P^ IHm iGIlTO IBlTZl |HLyZl^ we investigate the stability 

of profiles (u^,z £ ) in the ZND limit e — > for a class of models (jl.ip including both the 
Majda mode0 studied in p i [Lj |RV ] \LyZ2\ ILRTZ] and the physical (rNS) equations studied 



in P [GSl |Wl [ZH [LyZTj [JLWllTZi] . Our conclusion, confirming a conjecture of [LRTZj . is 
that (linear and nonlinear) stability in the ZND limit is equivalent to viscous stability of the 
component Neumann shock profile together with hyperbolic stability of the associated ZND 
detonation. 

Together with the results of |HLyZl| verifying viscous stability of ideal gas shocks, 
this gives a rigorous connection between viscous stability of (rNS) detonations and inviscid 
stability of the associated (ZND) detonations, yielding immediately a number of stability 
and bifurcation results through the extensive (ZND) literature. Specialized to the Majda 
model, it recovers the sole previous result, due to Roquejoffre and Vila [RVJ. 



1.1 Assumptions 

Loosely following [ZH E21 IMaZ31 IMaZ41 ILRTZllTZi] . we make the assumptions: 
(HO) /, B, cb, C G C 2 . 

(HI) The eigenvalues of df (u) are real, distinct, and different from s, for all u near 
the image of ZND profile u° , in particular for u = n_ , u* , u+ . 

1 Strictly speaking, a variant [L] with nonzero z-diffusion C > 0; however, the extension to the original 
Majda model is straightforward, substituting the weighted norm analysis of Sattinger [Sa] for the pointwise 
analysis of LRTZ] in order to conclude nonlinear stability. 
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(H2) B = f jj °\ with Slab, ftaC > 9 > 0, and d/ = f^j , d/n and d/i 2 

constant, with the eigenvalues of dfu real, semisimple, and of one sign relative to s, for all 
u under consideration (i.e., near a given detonation profile). 



(H3) ^ - i£df(u) - £ 2 B(u)j < 9>0, for all f G R, and for all u near the 

image of ZND profile u , in particular for it = u_, it*, u+. 

Remark 1.1. i?y 6/ocA; upper-triangular structure, we obtain from (H3) also 

(i.io) r*w °Vf 2 ( B i u) ° .w; q it\))<- ,e 



Oy V c (u,z)J \0 -k(j>{u)j ) ~ l + £ 2 ' 

> 0, /or a// £ G R, and /or a// u near the image of ZND profile u , in particular for 
u = u_, it*, an assumption in the nonlinear stability /bifurcation analysis of [TZJ^]. 

Regarding connecting profiles, we make the further assumptions: 

(PI) There exists a ZND profile u(x, t) = u°(x — st) of (|1.6p . smooth for x ^ 0, with 
u°(0~ = it*) and n°(0 + = it+), that is transversal in the sense that df(u°) — si is invertible 
for all x, so that the profile is locally unique by (|1.8p , 

(P2) There exists a viscous Neumann shock profile 

(x — st \ 
), lim u(x) = u*, lim u(x) = u + 

of the associated nonreacting Navier-Stokes equations ut + f{u) x = e(B(u)u x ) x , i.e., a 
connection between it*, u + of the traveling-wave ODE B(u)u' = f(u) — f(u+) — s(u — u+)± 
that is transversal in the sense that dfu(u) — s (constant by assumption (H2)) is invertibleo 

(P3) For 5 > fixed and e > sufficiently small, there exist viscous detonation 
profiles u e of (jl.ip . (|1.2p satisfying for some C, 9 > 0, and < k < 2, 

(1.12) |^((u e ,z e ) - (u°,z°))|(x) < Cee~ e W for x < -5, 

(1.13) |9f((n £ (x),z £ ) - (n(x/e),l))| < Ce + Ce 1 -^-^ 1 / 6 for - 5 < x < 0, 
and 

(1.14) |d£((u e (x),z e ) - (u(x/e),l))| < Ce 1 "^-^ 1 ^ for x > 0. 



2 This implies that the traveling wave ODE is nondegenerate type; the profile is then necessarily transver- 
sal by the extreme shock assumption (|l,5[l MaZ3 . 
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Example 1.1. The physical single- species reactive compressible Navier-Stokes equations, 
in Lagrangian coordinates, are \Ch\ \TZJ$ 



(1.15) 



d t r - d x u = 0, 
d t u + d x p = d ai (i/r~ 1 d x u), 
d t E + d x {pu) = d x (qdT~ 2 d x z + kt~ x 3 x T + vt~ 1 u8 x u). 
{d t z + k^(T)z = d x (dr- 2 d x z), 



where r > denotes specific volume, u velocity, E = e + ^u 2 + qz > total specific energy, 
e > specific internal energy, and < z < 1 mass fraction of the reactant. Here, v > 
is a viscosity coefficient, n > and d > are respectively coefficients of heat conduction 
and species diffusion, k > represents the rate of the reaction, and q is the heat release 
parameter, with q > corresponding to an exothermic reaction and q < to an endothermic 
reaction. Finally, T = T(r, e, z) > represents temperature and p = p(r, e, z) pressure. 

Under the standard assumptions of a reaction-independent ideal gas equation of state, 
p = Tr~ 1 e, T = c~ 1 e, where c > is the specific heat constant and T is the Gruneisen 
constant, and a smooth ignition function (j) vanishing identically for T < Ti and strictly 
positive for T > Ti, it is shown in JMaZ3[ \ TZJ$ that each of (H0)-(H3) are satisfied. 
Likewise, (Pl)-(PS) have been verified in IHuZB. in this case. 

Remark 1.2. We expect that (PS) can be shown by an argument like that of to be 
a general consequence of (H0)-(H3) and (P1)-(P2) and not an independent assumption. 
Note that (P3) typically requires more regularity than we assume in (HO). 

Remark 1.3. In (HI), it is enough that eigenvalues be semisimple. In the present, spectral, 
analysis, we use only that eigenvalues are real and distinct from s (used to separate decaying 
and growing slow modes in the analysis of Section \6.2.0) . In the linearized and nonlinear 
stability analysis for viscous detonations, one may relax the strict hyperbolicity assumption 
of \TZ~4J to semisimplicity, as discussed for the viscous shock case in \MaZ4 , IZl]/ . 



1.2 Main results 

Recall that, associated with the linearized eigenvalue problems for ZND and rNS detona- 
tions are the Evans-Lopatinski determinant Dznd and the Evans determinant D^, NS , each 
analytic on > — rj < 0, with zeros corresponding to normal modes of the respective 
linear problems; see Sections [2] and [3] for precise definitions. Likewise, there is an Evans de- 
terminant -Dtvs associated with the linearized eigenvalue problem for the associated viscous 
Neumann shock of the nonreacting Navier-Stokes equations (NS); see Section [6l 

Weak Evans-Lopatinski stability of ZND detonations is defined as nonvanishing of 
Dznd on KA > and strong Evans-Lopatinski stability as nonvanishing on 5RA > ex- 
cept for a simple zero at A = |ErH IEr2| \Z1\ IJLW| . Similarly, weak Evans stability of rNS 
detonations is defined as nonvanishing of E>^ NS on KA > and strong Evans stability as non- 
vanishing on IRA > except for a simple zero at A = [ZU \LyT\\ [L^Z2l iJLWl ILRTZllTZi] . 
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Likewise, weak Evans stability of the associated viscous Neumann shock is defined as non- 
vanishing of Dns on 3f?A > and strong stability as nonvanishing on 3f?A > except for a 
simple zero at A = 0: equivalently, nonvanishing of Djv f (A) on KA > |MaZ3L IZll IZ2l IZ3] . 

The following result established in [LRTZ1 ITZ4] equates strong Evans stability with 
linear and nonlinear stability of rNS detonations. A corresponding result holds for the 
component viscous Neumann shock |MaZ41 IZ2l El EEJ iHRZl [RZllZij . 

Proposition 1.4 ( |TZ4| ). Given (HO)-(HS), a viscous detonation profile (jl.2p of (II. ip is 

L 1 n LP — > LP linearly orbitally stable, p > 1, if and only if it is strongly Evans stable, in 
which case it is L 1 (1 — > L p PI H 3 asymptotically orbitally stable, for p > 1, with 

\(u, z)(-,t) - («, z){- -st- a(t))\ L v < C\(uq, zo) - (u, z)\ L i nH3 (l + t)~^ x ~p\ 
(1.16) \a(t)\ <C\U^-U £ \ L i nH3 , 

\a(t)\<C\US-U e \ L x nH *(l + t)-* 

for some a(-), where (u,z) is the solution of (jl.ip with initial data (uq,zq). 

Proof. This was established in Theorem 1.12, |TZ4j . for the case described in example ll.il 
However, the proof relied only on (H0)-(H3) and the consequent (ll.lOp . hence extends to 
the general case. □ 



Our main theorem is the following result linking Evans stability of (rNS) profiles (u e 
in the limit as e — > with Evans-Lopatinski stability of the limiting (ZND) profile (it , 



0\ 



Theorem 1.5. Assuming (H0)-(H3), (P1)-(P3), weak Evans stability of (u £ ,z £ ) for all 
e > sufficiently small implies weak Evans stability of the viscous profile of the compo- 
nent Neumann shock together with weak Lopatinski stability of the limiting ZND detonation 
(vP,z°), while strong Evans stability of (u e ,z e ) for all e > sufficiently small is implied by 
strong Evans stability of the viscous profile of the component Neumann shock together with 
strong Lopatinski stability of the limiting ZND detonation (u°,2p). 

More precisely, (i) For e > and r] > sufficiently small, there are no zeros of D r N$ on 
5RA > —rj for |A| > C/e, C sufficiently large, (ii) For C <|A| < C/e, C sufficiently large, 
on UiA > —77 for 77, e > sufficiently small, e times each zero of D^ NS converges to a zero 
of Dn s( x ^ on 3ftA > 0; moreover, each zero of D^s on IRA > is the limit of e times a zero 
ofDrNS on J?A > ; for C < |A| < C/e. (Hi) For |A| < Co, Co arbitrary, on KA > -77 < 0, 
the set of zeros of D^ NS converges as e — > to the set of zeros of Dznd, for any sufficiently 
small 77 > such that Dznd does not vanish for 3?A = —77 and |A| < Co. 

Proof. Assertions (i)— (iii) are established in Proposition 17.11 Corollary 16.21 and Corollary 
15.31 whence the remaining assertions follow by definition of weak and strong stability of the 
various waves. □ 

Remark 1.6. Assuming Evans stability of the associated viscous Neumann shock, we re- 
cover from (ii)-(iii), taking Cq — > 00 and e — > 0, the somewhat delicate result of \Erl\ \Er£tf 
that Dznd does not vanish on KA > for |A| sufficiently large. 
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Remark 1.7. Recall \JLW\j that, for fixed e, low-frequency Evans stability, defined as strong 
Evans stability for |A| < cq, some cq > 0, is equivalent for either ZND or rNS detonations to 
the simpler condition of "Chapman- J ouget" stability; see JJLW)/ for further details. Thus, 
low-frequency stability of ZND waves is necessary for Evans stability of rNS detonations, 
along with weak Evans-Lopatinski stability as stated in Theorem \1.5[ 



1.3 Discussion and open problems 

Together, Proposition 11.41 and Theorem 11.51 give a rigorous connection between Evans- 
Lopatinski stability of (ZND) detonations and nonlinear stability of nearby (rNS) detona- 
tions for e > sufficiently small, giving a satisfying mathematical validation of physical 
conclusions made through the extensive ZND stability studies in detonation literature. By 
contrast, to our knowledge there is no analog of Proposition 11.41 for the (ZND) equations 
themselves, and indeed the physical meaning of Evans-Lopatinski stability in that context 
in terms of nonlinear stability or well-posedness is unclear; see |JLW| for further discussion. 

More, convergence of zeros of D e rNS to those of Dznd implies that finer phenomena such 
as Hopf bifurcation are inherited in the ZND limit as well as stability. This is perhaps more 
important, as it is well-known that ZND detonations are frequently unstable, bifurcating 
to pulsating and cellular fronts. See [TZ2] ITZ3|. ITZ4|, ISSj IBeSZj for a rigorous discussion of 
such bifurcations in the context of (rNS). 

Existence, stability, and bifurcation of detonations away from the ZND limit are im- 
portant open problems. Such general situations appear to require numerical investigation, 
as is standard in the combustion literature even for the simpler (ZND) model; see, e.g., 
[LSI IHuZ2] and references therein. Treatment of the ZND limit in multi-dimensions is an- 
other important open problem. In particular, the implications for nearby (rNS) profiles of 
high-frequency ZND instabilities pointed out in [Er3j is an intriguing mathematical puzzle; 
see |JLW| for further discussion. 

We note that one-dimensional stability in the ZND limit was previously established in 
[RV| for detonation profiles of Majda's model. Our results both recover and illuminate 
this prior result, since the associated Neumann shock profile, because it is scalar, is in this 
case always stable (see, e.g., [Saj), as is the limiting ZND detonation. Related singular 
perturbation results for systems of conservation laws may be found in |PZl |HLZ, CHNZ, 
HLyZl, HLvZ2 IBHZ| and (for multi-wave patterns) in [QZ1[ [Z~3] . In particular, we rely 



heavily on the basic methods of analysis developed in [PZj, [HLZj , |Z3j . and [BHZJ. 



2 The Evans-Lopatinski determinant for (ZND) 

We begin by recalling the linearized stability theory for ZND detonations following [Erl, 
IZlt IJLWj IHuZ2] . Shifting to coordinates x = x — st moving with the background Neumann 
shock, write (jl.6p as 

(2.1) W t + F(W) X =R(W), 
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where 

<»> »"-(:). «-(!t:t) 

To investigate solutions in the vicinity of a discontinuous detonation profile, we postulate 
existence of a single shock discontinuity at location X(t), and reduce to a fixed-boundary 
problem by the change of variables x — > x — X(t). In the new coordinates, the problem 
becomes 

(2.3) W t + (F(W)-X / (t)W) x = R(W), x^O, 
with jump condition 

(2.4) X'(t)[W)-[F(W)) = 0, 

[h(x,t)} := h(0 + ,t) — h(0~,t) as usual denoting jump across the discontinuity at x = 0. 



2.1 Linearization 

In moving coordinates, is a standing detonation, hence (W°,X) = (W°,0) is a steady 

solution of |[%3]H[%3|) . Linearizing (pO ]) - (pl4"]) about (W"°,0), we obtain the linearized 
equations 

(2.5) (W t - X'(t)(W°)'(x)) + (Aiy)j = EW, 

(2.6) I'(t)[f°]- W=0, x = 0, 
where A := {d/dW)F, E := (d/dW)R. 



2.2 Reduction to homogeneous form 

As pointed out in [JLW ] . it is convenient for the stability analysis to eliminate the front 
from the interior equation (|2,5p . Therefore, we reverse the original transformation to linear 
order by the change of dependent variables 

(2.7) W - X(t)(W°)'(x), 

motivated by the calculation W(x - X(t),t)) - W(x,t) ~ X(t)W x (x,t) ~ X(t)(W°)'(x), 
approximating to linear order the original, nonlinear transformation. Substituting (|2,7p 
in (|2.5p - (|2.6|) . and noting that x-differentiation of the steady profile equation F(W°) X = 
R(W°) gives (A(W )(W )'(x)) x = E(W )(W°Y (x), we obtain modified, homogeneous in- 
terior equations 

(2.8) W t + (AW) X = EW 

agreeing with those that would be obtained by a naive calculation without consideration of 
the front, together with the modified jump condition 

(2.9) X'(t)[W°) - [A(W + X{t)(W )')} = 
correctly accounting for front dynamics. 
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2.3 The stability determinant 

Seeking normal mode solutions W(x,t) = e xt W(x), X(t) = e xt X, W bounded, of the 
linearized equations (|2.8p - (|2.9p . we are led to the generalized eigenvalue equations 

(2.10) (AW)' = (-XI + E)W, x^O, 

X(X[W°] - [A(W )']) - [AW] = 0, 

where "/" denotes d/dx. or, setting Z := AW, to 

(2.11) Z' = GZ, x + 0, 

(2.12) X(X[W°] - [A(W )'}) - [Z] = 0, 
with 

(2.13) G := (-XI + E)A-\ 

where we are implicitly using the fact that A is invertible, by (PI) and s > 0. 

Lemma 2.1 ( (ErTj [Er2 l [JEW] ) . On RX > 0, the limiting (n+ 1) x(n+l) coefficient matrices 
G± := lim z _ > ± o G(z) have unstable subspaces of fixed rank: full rank n + 1 for G+ and rank 
n for G-. Moreover, these subspaces extend analytically to RX < — rj < 0. 

Proof. Straightforward calculation using upper-triangular form of G± \Erl\ IEr2| \Z1\ iJLWj . 

□ 

Corollary 2.2 (jZT jlJLWj ). On RX > 0, the only bounded solution of (j2~TTj) for x > is 

the trivial solution W = 0. For x < 0, the bounded solutions consist of an (n) -dimensional 
manifold SpanjZ-^ , . . . , Z+}(X, x) of exponentially decaying solutions, analytic in X and 
tangent as x — > — oo to the subspace of exponentially decaying solutions of the limiting, 
constant- coefficient equations Z' = G~Z; moreover, this manifold extends analytically to 
RX < -T] < 0. 

Proof. The first observation is immediate, using the fact that G is constant for x > 0. 
The second follows from standard asymptotic ODE theory, using the conjugation lemma of 
[MeZl| (see Lemma |A.1| Appendix [A]) together with the fact that G decays exponentially 
to its end state as x — > — oo. □ 

Definition 2.3. We define the Evans-Iopatinski determinant 

D ZND (X):=det(Z^(X,0), Z~(X,0), X[W°] - [A(W )']) 

= det(Zf(A,0), Z-(X,0), X[W°] + A(W )' (()-)) , 



where Z- (X,x) are as in Corollary ] 2. E 

The function Dznd is exactly the stability function derived in a different form by Er- 
penbeck |ErTJ IEr2] . The formulation (|2.14p is of the standard form arising in the simpler 
context of (nonreactive) shock stability |Er4] . Evidently (by (|2.12p combined with Corollary 
12. 2[) . A is a generalized eigenvalue/normal mode for MA > if and only if Dznd(X) = 0. 



10 



3 The Evans determinant for (rNS) 
3.1 Linearization 

Linearizing (|l.ip about (u e , z £ ) in moving coordinates yields linearized eigenvalue equations 

(3.1) \W + {AW) X = EW + {eBW x ) x , 

^=(«) = (^.^ = (:).*- = (?). s = (Ss 

and 

(3.2) Av := dF{W e )v - e(dB(W e )v)W£, (An, An) = (df n ~ s, (df 12 ,0)) = constant, 
E, F as in OT8|) . 




3.2 Expression as a first-order system 

Setting x = x/e and 



Y \ (AW - eBW x \ (Y{\ „, ^u 2 



we may write (J3JJ) as a first-order system 

(3.3) W = g e (X,x)W, 
where 

(e{E - A)i n X -e(£ - \)A£A 12 

(3.4) £ £ := e(E — A) 

VB^AaiAJi 1 -H 6- 1 (i 2 2-i 2 iiri 1 ^i2)y 

and " denotes d/dx. Here, we are using implicitly the facts that An and b are invertible, by 
CT - (P2), and (H2). 

We have the following analogs of Lemma l2.1l and Corollary l2.2l which follow in essentially 
the same way; see fLRTZl ITZ4| for further details. 

Lemma 3.1 ((ZTJ ITZ4| ). On RX > 0, the limiting (n + 1 + r) x (n + 1 + r) coefficient 
matrices Q± := lim^-too ^(z), r = dimu2 + 1, /tawe stable subspaces of fixed, equal rank r. 
Moreover, these subspaces extend analytically to IRA < —n < 0. 

Corollary 3.2 ((ZT j [TZ4] ). On MA > 0, the bounded solutions of flU]) on x < consist 
o/ an (n + 1)- dimensional manifold Span{W-f , . . . , W~, 1 }(A,x) o/ exponentially decaying 
solutions in the backward x direction, analytic in A and tangent as x — > — oo to the subspace 
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of exponentially decaying solutions in the backward x direction of the limiting, constant- 
coefficient equations W = Q-W; the bounded solutions of (|3.ip on x > consist of an 
(r — 1)- dimensional manifold Span{>V„ + 2 , . . . , Wi_ 1+r }(A, x) of exponentially decaying so- 
lutions, analytic in A and tangent as x — > +00 to the subspace of exponentially decaying 
solutions of the limiting, constant- coefficient equations W = G+W; Moreover, these mani- 
folds extend analytically to IRA < — 77 < 0. 



3.3 The stability determinant 
Definition 3.3. We define the Evans function 

(3.5) D^ 5 (A):=det(V\^, W~ +1 , W+ +2 , W+ +1+r ) (A, 0), 

where W- (A, x) are as in Corollary \3.2l 



4 Fast vs. slow coordinates 

Note that the coordinate transformation to stretched, or "fast" variables 
(4.1) (x, t, A) = (x/e,t/e, eX), 

changes equations (jl.ip to 

^ 4 ^ U't + fi u )x = {B(u)uz)x + ekqip(u)z, 

z- t = (C(u, z)zx)s ~ ekip(u)z, 

shifting the small parameter e from diffusion to reaction terms. The computations to follow 
may be thought of as alternating between (original) slow variables and fast variables in our 
analysis, as convenient for different regions in A and x. In particular, first-order equations 
(|3.3p ~ (|3.4p may be recognized as the first-order linearization of fast equations (|4.2|) . 



5 Region I: |A| < C 



We first study the critical "ZND" region |A| < C, or e/C < |A| < Ce, where behavior of 
the (rNS) Evans function D r ^s is governed by that of the Evans-Lopatinski determinant 
Dznd- Setting M » 1 to be a large constant to be determined later, we study separately 
the zones x = x/e < —M and x = x/e > —M, on which the profile W £ is dominated 
respectively by the (ZND) profile W° and the viscous shock profile W, as described in (P3). 

5.1 "Slow", or "reaction" zone, x < —M 

Note that A, hence N := b~ 1 (A 2 2 - A 2 iA^A 12 ) is invertible for x < -5, by ([33]) and 
(P1)-(P2). Setting W = TZ, where T := (_^-i £ }) and i := (ft-^l^nS 

-N-H = -(A 22 - A 21 A^A 1 2)- 1 (A 21 A^, -I), 
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and noting that N, I ~ W £ on x < -M, by (P3), we transform ([Oft - lpTl]) to 

(5.1) z = n £ z, 

(5 2) H e - T~ l G £ T - T~ l f - ( e{E ~ X ^ £m \ + ( . °\ 

(5.2) W ^ iV + 0(£)J + ^C>( £ + |W~|) oj' 

where m := {E — A) ( ^ J, and terms O(-) are smooth and analytic in A. Here, 

N, by (|1.5p and the relation between viscous and inviscid shock structure [MaZ3| , has one 
positive eigenvalue and r — 1 negative eigenvalues, where r = dimu 2 + 1. 

/j — gjV -1 m\ 

By a further transformation Z = SV, S = I , we transform to V = IC £ V, 







where 



(5.3) /C E = S-'H'S-S-'S = [ v v ^) + 



/0(e 2 +e|W^ e |) 0(e 2 + e|# e |; 
y 0{e + \W £ \) 0{e+ \W £ \) 



By a further transformation f ^ ^ if necessary, we may take without loss of generality 

N block-diagonal, AT = ( 1 V with KiVi > ry > and 5?N 2 < —77 < 0, both with a 

V u iv 2 y 

uniform spectral gap from block e{E — X)A~ l ~ e. 

Applying the asymptotic ODE results of Lemma I A. 71 with r/~l, 5 = 0(e + \W £ \), 
and |G| < C, and of Remark ESl with r) ~ 1, 5 := e(e + |# e |), |9i 2 | + |©u| < C, 

/ Wi 

I ©21 1 + I ©22 1 — C/e, we find that there is a change of coordinates Z = SU, U = fU^' 
defined on x < — M, 

(5.4) S- 1 =[ _ r = - 1+1 U = ol / ol 

close to the identity converting (|5.ip into three decoupled equations 

ii = e(£ - X)A~ 1 Zi + 0(e 2 + e|# £ |)^i; U 2 = 0,U 3 = 0, 

(5.5) Z 2 = (JVi + o(l))Z 2 - U x = 0,U 3 = 0, 
Z 3 = (N 2 + o(l))Z 3 ; U X = Q,U 2 = Q. 

Focusing on the "slow" Z\ mode, changing coordinates from x back to x = ex, we obtain 

(5.6) Z[ = (E- X)A- 1 Z 1 + 0(e + [W^x/e^Zt, 
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where "'" denotes d/dx, and all coefficients converge uniformly as Ce~ v ^, n > to their 
limits as x — > — oo. Applying the convergence lemma, Lemma IA.4I (Appendix IA.2|) . for 
x < —M, together with Remark I A. 51 (|A. 16|) . we find that there is a further coordinate 
change Z\ = PZ , P = I+0(s+e~ vM ) = I+o(l) for \x\ < 2M, taking the decaying/growing 
modes of (|5.6|) to those of 

Z' = (E — X)A- 1 Z, 

which may be recognized as exactly the interior equation (|2.1ip . (l2.13j) associated with the 
Evans-Lopatinski development. 

Tracing back through our coordinate transformations, we find that the n slowly-decaying 
modes W-f , . . . W~ as x — > — oo are given at x = —M by 

(5.7) W-=(I + o(l))( Z A, j = l,...,n, 

where Z~ are as in (I2,14p . while the single fast-decaying mode as x — > — oo is given at 
x = —M by 



(5.8) W~ +1 = c(-M) 



o(l)ev 
(J + (1)W' 



where c(x) is exponentially decaying as x — > — oo and v lies in the unique stable eigendirec- 
tion of N(—M), N := b~ 1 (A22 — AnA^Ax?). Here, o(l) depends on M, going to zero as 
M — * oo and e — > at the same time. 

Remark 5.1. At e = 0, the single fast-decaying mode reduces to the translational zero 
eigenfunction of the associated viscous Neumann shock (necessarily fast- decaying), and so 
we may refine (|5.8p to 

(5-9) W- +1 =( ° il)£ KY 



z 



5.2 "Fast", or "Neumann shock" zone, x > —M 

Applying the convergence lemma, Lemma [A.41 for x > —M, using the asymptotics of (P3), 
together with Remark IA.5|, (|A.16h , we find that there is a change of coordinates PI with 
P+ = l + 0(e) for \x\ < M, such that W := P+V converts (I3.3p - (l3.4p to the same equations 
with e = 0, i.e., the shock eigenvalue system at A = 0, hence, by inspection, the decay- 
ing modes V^" +2 , • • • , V^~ +1+r at +oo, necessarily fast-decaying, by the Lax characteristic 
assumption (II. 5|) (see [Ma Z3] ) , are by inspection (see [Zl] for similar calculations) of form 





(5 - 10) V ^ = \vM 
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Vj independent, hence for \x\ < M 



(o.ll.) W+ = (/ + 0(e)) ( ° ) . ./' - " + 2. // + I + r. 



Here the constant O(-) depends on (growing exponentially with) the fixed constant M, 
since we are shifting x — > x + M in order to apply the convergence lemma. Evaluating at 
at x = —M, we have 



(5.12) W+(-M) = (/ + 0(e)) ( Vj{ °_ M) ), j = n + 2,...,n+l + r. 



5.2.1 Variation in e 



At this point, gathering information, and noting, by Abel's formula, that the Wronskian 
D r NsW evaluated at x = is equal to a nonzero constant times the same Wronskian 
evaluated at the point x = —M at which we have information about solutions from both 
sides, we have that D r jsrs(\) is proportional by a nonzero constant, analytic in A, e, to 

det ((/ + o(l)) (^) ,...,(/ + (1)) ( Z ; ) , (/ + o(l)(e) ( J + J , 

(5 - 13) (I + O(e))(„ ),...,(I + O(e))( ° 

= 0(e), 



V n +2 J \V n +l+r 



since up to 0(e) there are only n nonzero entries Z- in the (n + l)-dimensional Y block. 
Here, the constants 0(e) and o(l) by construction are analytic in A, e as well. 

fd~u\ 

This reflects the fact that at e = there is a solution W = dxW = [ x I decaying at 



v / 

both ±oo of (|3.ip . corresponding to the translational eigenmode of the linearized equations 
about the associated viscous Neumann shock; see Remark 15.11 To extract the next-order 
behavior, we compute the first variation of this special mode with respect to e at e = 0, by 
an argument similar to those used in \G7\ IZS} |LyZl[ |LyZ2| to compute the first variation 
with respect to A at A = 0. 

Specifically, recalling that u n +l> • • • , v n+ i +r are a basis in C r , we may perform a column 
operation using the final r columns to cancel the entry v n +i in the (n + 2)nd column in 
determinant (15.130 . to obtain 

(5.14) D £ rNS (\)=^(X,e)(I + o(l))edet(v n+2 ... v n+1+r ) det (Z~ ... Z~ -Y £ ) 
hence 

(5.15) = det (Zf ... Z- -Y £ )+o(l), 
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as e — > 0, o(l) — > as M — > 00, where ip and ^ := ipdet (v n+ 2 ■ ■ ■ v n +i +r ) are nonvan- 
ishing factors analytic in A, e, and 

(5.16) Y £ :=d £ Y*\ £=0 =d £ Wt\ £ =o, 

where W*(e,A, x) is the (necessarily fast-) decaying solution of (|3.ip at x — > +00 defined 
by 

(5.17) W = P£(A,x)( ° 
where Pl(X,x) is the conjugating transformation described above, and 

W *\ £=0 = d -W=( d f 

Writing (|3.ip in x coordinates as 

(5.18) e(E - X)W = (A £ W - B £ Wx)~ x = Y, 

and differentiating (|3.ip with respect to e, we obtain the variational equations 

(£ - \)W*\ £= o - (d £ A £ W* - d £ B £ W£) £ \ £=0 = Y £ , 

or 

(5.19) EW S - (XW + d £ A £ W~ x - d £ B £ W- x - x ) x = Y £ . 
Finally, recall that differentiating the traveling-wave ODE 

(AW £ ) X = eEW £ + (BWi) s 

with respect to x yields 

(5.20) EW% = e-\AW £ x - BW,,, ;,. = (AW £ - eBW £ x ) x . 
Together with the estimates 

(5.21) |W|- W s \ < Cee-iW, 
r] > 0, for x > —M and 

(5.22) \Wi(-M)\, \W XX (-M)\ =o(l) 

coming from asymptotics (P2), (P3), we find, combining (|5.19p . (|5.20p . (|5.2ip . and (|5.22p 
that 

(Y £ ) x = E{W X - W§) + (AW £ - eBW £ xx ) x - (XW + d £ A £ W x - d £ B £ W xx ) x 
' ' ! '"" ' = (AW £ — XW — sBW £ x - 8 £ A £ W X + d £ B £ W xx ) x + Oiee'^). 
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Integrating (|5.23j) in x from x = —M to x = +00 thus yields 

(5.24) - Y e {-M) = {-AW e x + XW + eBW £ xx + d e A e W x - d £ B e W xx )\l= + _^ + O(e), 
From asymptotics (P3), we find that 

(5.25) AW% - eBW e xx = AW0 X + (AW X - eBW xx ) + o(l). 

Recalling (by differentiation in x of the traveling- wave ODE for the viscous Neumann shock 
profile) that 

AW X - eBW xx = 0, 

we find, substituting (|5.25p into (j5.24jl and discarding lower-order terms and vanishing 
boundary terms at +00, 

(5.26) - Y £ {-M) = X[W°] - [A(W°)'{0)} + o(l), 

where W° as in (PI) denotes the associated ZND profile and [■] the jump in values across 
its Neumann shock discontinuity. 

5.3 Convergence to D Z nd 

Proposition 5.2. Assuming (HO)-(HS), (Pl)-(PS), for |A| < C and > -77, rj > 
sufficiently small, converges uniformly as e — > to Dznd{X), where ■) is a 

nonvanishing factor that is analytic in A. 

Proof. Choosing monotone sequences Mj (increasing) and £j (decreasing) such that o(l) < 
l/j in (|5.26p and (|5.15p . define §(A,e) to be equal to the function *(A,e) in (|5.15p that is 
associated with Mj, where j is the maximum integer such that e < Ej. Then, ^ is analytic 
in A by construction, and, combining (|5.15p . (|5.26p . and the definition of Ej, Mj, we have 



DzndW <C/j^0 ase^O. 

e*(A,e) 



□ 



Corollary 5.3. Assuming (HO)-(HS), (Pl)-(PS), for |A| < C and lk\ > -rj, r] > 0, the 
set of zeros of D e rNS converges as e — > to the set of zeros of D^ NS , for any sufficiently 
small rj > such that Dznd does not vanish for 9?A = —rj and |A| < C . 

Proof. Noting that zeros of Df. NS agree with zeros of we obtain the result by 

properties of uniform limits of analytic functions. □ 
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6 Region II: C/e > |A| > C » 1 

We next consider the "Neumann shock" region Cs < X < C, C » 1, in which behavior of 
D t ns is dominated by that of the Evans function Dj^g of the associated viscous Neumann 
shock. Here, D^s is defined similarly as in Definition 13.31 of D v ns as 

(6.1) £fr(A):=det(Wr, VV~ +1 , W+ +2 , W+ +1+r ) (A, 0), 

where VV^ are decaying modes at ±oo of the linearized eigenvalue equation 

(6-2) XW + {AW)x = {BWx)x 

about W, written as a first order system 

-Ai^i 1 -XA^A 12 
(6.3) W = g(X,x)W, G-\ —A 



K b- 1 A 21 A^ -b- 1 b- 1 (A 22 -A 21 A^A 



12) 



withW= (X) = \ AW - BW ) , B := B(W), and Av := dF{W)v-(dB(W)v)W x , where 



F is as in (|2.8p . For further details, see, e.g., [ZH IZ2] 



6.1 Fast zone x > —M 

Noting that A = eA is bounded by assumption, we have by (|3.4|) and (P3) 

/eE-XA^ -eE - ~XA^A 12 \ 

(6.4) G e =\ eE-X \=g + 0(ee^ x \), 

\b- 1 A 21 A^ 1 1 -b- 1 l' x {A 22 - A 21 A^A l2 )J 

where G is bounded and converges at uniform exponential rate to its limits at =Loo. Applying 
the convergence lemma, Lemma \AA\ for x > —M, together with Remark I A. 51 (|A.14|) . 
similarly as in Section 15.21 but now treating A as a fixed parameter, we find that there 
is a change of coordinates Pf with Pf = I + 0(e) for x > —M, such that W := PfVV 
converts ()3.3p - (|3.4p to the viscous shock system ()6.3p . Evaluating at x = —M, we thus 
have W+(A, —M) = (I + 0(e))W+(A, — M), or, by the assumption that |A| >> e, 

(6.5) Wf(\,-M) = (l + o(X)) Wf(X,-M), j = n + 2,...,n + l + r. 

6.2 Slow zone x < —M 

6.2.1 Case a. C/e > \X\ > l/Ce, C > arbitrary 

We first treat the easier case l/Ce < |A[ < C/e, or C _1 < |A| < C, for arbitrary C > 0. 
From (H3), it follows that G(A, — oo) has no pure imaginary eigenvalues for 5RA > and 
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A ^ 0. By continuity, the same holds for Q(X,x) for x < —M, > and A bounded and 
bounded away from zero. 

By the assumption \jC < |A| < C, therefore, the stable and unstable subspaces of Q 
have a uniform spectral gap for all x < —M, hence, by standard matrix perturbation theory 
[K l \ZH\ [Z5 | RjMWZ5j . there exist smooth transformations T(Q) such that 

T~ X QT = , RMx > V > 0, mi 2 < -r] < 0. 

Making the change of variables W = TZ, we thus obtain 

Z={T- 1 GT-T- 1 f)Z=[(^l l ^+o(l))z 

for x < —M. Applying the tracking lemma, Lemma \A.7\ we find that the manifold of 
solutions of (I3.3p - ()3.4p decaying at — oo is within angle o(l) of the unstable subspace of 
Q(X,x) at each x < —M, in particular for x = —M. 

But, by the same reasoning, the manifold of decaying solutions of (16, 3D at x = —M 
is also within angle o(l) of the unstable subspace of Q, which, by continuity in e/uniform 
spectral gap is within angle 0(e) = o(l) = o(A) of the unstable subspace of Q. From this, 
and (|6,5p . it follows that, up to a normalizing factor ^>(e, A) that may be taken analytic in 
A, 

(6.6) -£^ = D Ns a) + o(X) 

on 5RA > —7], C _1 < |A| < C, uniformly as M — > oo, and e — > 0, for C > arbitrary. 
6.2.2 Case b. 1/Ce > |A| > C » 1 

Finally, we treat the more delicate case C < |A[ < 1/Ce, or Ce < |A| < C _1 , with C >> 1. 
Proceeding as in Section 15.11 by the series of coordinate transformations (|5.ip - (|5.6p . but 
taking account of the different order of A in this case, in particular, noting that em in (15. 2j) 
is now O(A) and not 0(e) as before, we obtain 



,0.7) ]C e = ((eE-X)A-i + X>(3 ) + 



/o(e|A| + [A[|# £ |) 0(e|A| + [A||# £ |)\ 

y 0(e+|# e |) 0(e+[# e |) y 

in place of (15.3H . where /? is the "frozen-coefficient corrector" obtained by dropping terms 
involving x-derivatives of the transformations involved, 





o(\x\) o(\x\y 

I o(l 
o(l) I 
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for x > — M in place of (|5.4p . and 

(6.9) ^ = 

in place of (|5,8p . where u lies in the unique stable eigendirection of N{—M), with 

iV:= b-^iaa-iaiA^iia). 
By the same reasoning, applied to the viscous shock equations (|6.3p . we have 

o(l)\\\v 



W n+1 = c(-M) 



(I + o(l))v 



where v lies in the same unique eigendirection. and VV n +i is the single fast-decaying mode 
of (|6.3p at — oo. Comparing, we thus have, for some normalizing factor ip(M, e, A), analytic 
in A, 

(6.10) W" +1 (-Af) = [l + o(\\\))i>(M,e,\)W- +1 (-M). 

We now turn to the description of the remaining, slow-decaying, modes Wf,... W~. 
Focusing on the region x > — Ci\ log |A||, Ci » 1, we find that the slow, first, equation of 
(15.51) becomes 



(eE - A)i" x + f3X 2 jZ 1 + 0(e + |W £ (x/e)|)|A|Zi 
(eE - X)A~ 1 + f\ 2 )z 1 + o(l)|A| 2 Zi, 



(6.11) 



where ((sE — A)^^ 1 + f3X ) agrees to second order in A << 1 with the restriction of Q to 
its slow subspace, which, recall, is the direct sum of n slow-decaying modes as x — > — oo 
and a single slow-growing mode as x — ► — oo. Here, o(l) — > as C — > oo and e — > 0. 
Computing explicitly, and noting that z ~ e _ee ' x ', we have 



eE = e 



f qkd(p(u £ )z qk(f)(u) \ _ [0 qk<p(u)\ ~ 
\-kd<p{u e )z -k<f>(u)J \0 -k(f>(u)J ^ a[X,Ah 



where a = 0(ee~~ 9e ^) is both o(A) and uniformly bounded in L l (x). Thus, up to 0(a) + 
o(|A| 2 ) the eigenvalues of ((eE — A)A _1 +/3A 2 ) agree with the growth rates fj,j of slow modes 
e H x Wj, Wj = constant, of the nonreacting frozen-coefficient operator 

«-> »+rr' c( „°>i+(o q t$, 

obtained by neglecting 0(dxW £ ) derivative terms and setting z to zero. 

A standard low-frequency matrix perturbation computation [Zl} LyZl[ |LyZ2[ ITZ4|. 
ILRTZ] shows that the modes Wj of (|6.12p are analytic in A and e, lying approximately 
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in the eigendirections of A, with the associated growing (i.e., negative real part) eigen- 
value Hn+i separated in modulus by order |A| from the decaying onesjl moreover, by (jl.lOp . 
the growing eigenvalue is separated in real part by order ~ |Aj 2 from associated decaying 
(positive real part) eigenvalues. 

From this, and the fact that the entire coefficient matrix is order A, it follows that there 
exists a smooth matrix- valued function Q(W £ ) with 

Q->( (£ £-A)i- + m ; +a2 ), 

(6.13) »mi > f)\\\ 2 > 0, Km 2 < -r?|A| 2 < 0, (ayU^) < C, fj = constant. 
Making the change of coordinates Z\ = Qz, and noting that Q' = 0(W £ ), we thus obtain 

z= {Q-\{eE-~\)A- 1 +(3\ 2 )Q-Q- l Q)z= f™ 1 ^z+{o(\\\ 2 )+0(W £ ))z, 

with gap condition f|6. 131) . 

Since \W £ \ = 0(e) = o(A), and rrij are spectrally separated by modulus ~ |A|, there is 
a further smooth coordinate change z = TZy with TZ = I + o(l) converting the equations to 

dy = ctj + 0(W £ ), (3 = o(W £ )), where \ctj\ L i < C, \f3\ L i^)\ = o(l), and rrij satisfy (|6.13p . 
That is, we have an equation of form (|A.27p with 5 = o(|A| 2 ) and 7] ~ |A| 2 , so that S/rj = o(l). 

Applying the tracking lemma, Lemma IA.71 as generalized in (1A.28|) , Remark IA.101 at 
x = — C | log | A 1 1 , and untangling coordinate changes, we thus find that the slow-decaying 
modes W~, j = 1, . . . , n lie within angle o(l) of the stable subspace of Q{— C\ log |A| |), which 

in turn lies within angle o(l) of the stable subspace of G(A, — oo) and (by a repetition of 
the same argument), analytic multiples of the slow-decaying modes VV~, j = 1, . . . ,n, at 
x = -M. 

Finally, going back to the original equation (|6.11|) . and noting that Z\ = 0{\\\)Z\, we 
find that the change in Z\ in evolving from x = — C \ log |A|| to x = — M is order 

(e |A|| icslAM _ i;) |^ ri( _ C7[ i g | A ||)[ ~ |A|| log IAIM^xC— C| log |A||)| = o(l)|^!(— C| log |A||)|, 

hence at X = —M also the slow-decaying modes W~, j = 1, . . . ,n lie within angle o(l) of 
analytic multiples of the slow-decaying modes VV7 at x, = —M. Collecting facts, we have 

(6.14) Wr(-M) = {I + o(l)W j (M,C 2 ,e,\)Wr(-M), j = l,...,n, 

3 This does not require strict hyperbolicity of df, but only det(d/ — s) 0; see Remark 11.31 
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where ipj are nonvanishing and analytic in A. (Here, we are using also the fact that, by 
(|6.7p ~ (|6.8p . the manifold of all slow modes, both growing and decaying, stays angle o(l) 
close to the slow subspace of Q(X,x) for all x > —M.) 

Finally, collecting estimates (|6.5p , (|6.10p , and (|6,14p , we see that fast modes VV„- , up to 
nonvanishing analytic factor ipj, are given by (/+o(A)) times the correspong fast modes Wj 
of the associated shock stability problem, while slow modes W7 are are given by (/ + o(l)) 
times the correspong slow modes W~ . Recalling, similarly as in estimate (|5.13p . that at 
A = 0, the fast mode is a linear combination of the fast modes Wj, so that at A 

it remains within angle O(A) of their combination, we find, applying a column operation 
cancelling W n+1 to order A and factoring out A from that column, that we reduce the (rNS) 
determinant (]3.5p to 

(6.15) \(l + o(l))(U J $ j )D N s(X). 

The key observation here is that because only fast modes are involved in the vanishing of 
Dns at A = 0, only fast modes must be estimated to the sharper relative error o(A) in order 
to obtain the result (|6.15p . with o(l) tolerance sufficing for slow modes. 

Thus, we obtain again an estimate rTTTT^v; = ^nsW + on ^ > ~~ Vj Ce < 

|A| < C^ 1 , where o(l) — > uniformly as M — > oo, C — > oo, and e —> 0, with := Hjipj 
nonvanishing and analytic in A. Note that we are using here the assumption that C >> 1, 
which was not needed in case a. 



6.3 Convergence to D N $ 

Proposition 6.1. Assuming (HO)-(HS), (Pl)-(PS), For C/e > |A| > C » 1, KA > 0, 

eA^(A^e) conver 9 es uniformly as s — ► to Dn ^( £X ) ? where ^(-, ■) is a nonvanishing factor 
that is analytic in A. 

Proof. From (|6.6p . (|6.15p . we have in either case a or b that 

(6.16) DrNsW =D NS (eX) + o(l)eX 

*(M, C,e, A) 

on 5?A > — rj, uniformly as e — > 0. where o(l) — > uniformly as M — > oo, C — > oo, 
and e — > 0. Choosing monotone increasing sequences Cj, Mj and a monotone decreasing 
sequence £j such that o(l) < l/j in (|6.16p . define ^(A, e) to be equal to the function ^(A, e) 
in (|6.16p that is associated with Cj , Mj , where j is the maximum integer such that e < e, . 
Then, ^ is analytic in A by construction, and, combining (I6.16P with the definition of £j, 

Mj, Cj, we have §f^y - £li §r^ < C/j ^ as e ^ 0, giving the result. □ 

Corollary 6.2. Assuming (HO)-(HS), (Pl)-(PS), for C < |A| < C/e, C sufficiently large, 
on 5RA > —r\ for r\, e > sufficiently small, e times each zero of D e rNS converges to a zero 
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of N ^ on ^ > 0; moreover, each zero of D^g on 9?A > is the limit of e times a zero 
ofD rNS on KA > 0, for C < |A| < C/e. 

Proof. Noting that zeros of D e rNS agree with zeros of ~^fp~~y> we obtain the result by 
Proposition 16.11 and properties of uniform limits of analytic functions. □ 

7 Region III: |A| > C/e, C » 1 

Finally, we consider the straightforward "hyperbolic-parabolic" region |A| > C/e, C >> 1, 
or |A| >> 1, on which zeros of D r ^s are prohibited stable by basic hyperbolic-parabolic 
structure/well-posedness of the underlying problem (jl.ip . 

Proposition 7.1. D r ^s does not vanish for |A| > C/e, C » I, 5RA > 0. 

Proof. In fast coordinates x, |A| >> 1, this follows by the same high-frequency analysis 
used in [MaZ3j to treat the viscous shock case, based on the tracking/reduction lemma, 
Lemma lA.71 See Propositon 5.2, |MaZ3] . or Proposition 4.33, [Z2] , □ 



A Asymptotic ODE theory 

A.l The conjugation lemma 

Consider a general first-order system 

(A.l) W' = AP(x,\)W 

with asymptotic limits A± as x — > =Loo, where p £ M m denote model parameters. 
Lemma A.l ( |MeZlj |PZ]). Suppose for fixed 9 > and C > that 
(A.2) \A p -A p ±\(x,\) < Ce- 6 ^ 

for x ^ uniformly for (A,p) in a neighborhood of (Xo), po and that A varies analytically 
in A and smoothly (resp. continuously) in p as a function into L°°(x). Then, there exist 
in a neighborhood of (Ao,Po) invertible linear transformations P+(x, A) = I + 0+(x, A) and 
P-{x, A) = I + Q p _(x, A) defined on x > and a; < 0, respectively, analytic in A and smooth 
(resp. continuous) in p as functions into L°°[0, ±oo), such that 

(A.3) \&l | < Cie" 9 " 1 * 1 for z ^ 0, 

/or any < 9 < 9, some C\ = C\{9,9) > 0, and the change of coordinates W =: P±Z 
reduces (jA.ip to i/ie constant- coefficient limiting systems 

(A.4) Z' = A P ±Z for x ^ 0. 
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Proof. The conjugators P± are constructed by a fixed point argument [MeZlj as the solution 
of an integral equation corresponding to the homological equation 

(A.5) P' = A p P - A P ± P. 

The exponential decay (|A.2p is needed to make the integral equation contractive with respect 
to L°° [M, +00) for M sufficiently large. Continuity of P± with respect to p (resp. analyticity 
with respect to A) then follow by continuous (resp. analytic) dependence on parameters 
of fixed point solutions. Here, we are using also the fact that (|A.2|) plus continuity of A p 
from p — > L°° together imply continuity of e^ x '(A p — A p ±) from p into L°°[0, ±00) for any 
< 6 < 6, in order to obtain the needed continuity from p — > L°° of the fixed point 
mapping. See also [PZ1IGMWZ5] . □ 

Remark A. 2. In the case that A is block diagonal or triangular, the conjugators P± may 
evidently be taken block diagonal or triangular as well, by carrying out the same fixed-point 
argument on the invariant subspace of (1A.5I) consisting of matrices with this special form. 
This can be of use in problems with multiple scales; see, for example, Thm 1.16, [BHZf. 

Definition A. 3 (Abstract Evans function). Suppose that on the interior of a set in 
X, p, the dimensions of the stable and unstable subspaces of A±(\) remain constant, and 
agree at ±00 ( "consistent splitting" \AGJ$ ). and that these subspaces have analytic bases 
Rj extending continuously to boundary points of fl. Then, the Evans function is defined 
on £1 as 

(A.6) D p {\) := det(P + R+,P + R+,p-R^,p-R 2 )\ x=0 , 

where P± are as in Lemma \A.1[ 

A. 2 The convergence lemma 

Consider a family of first-order equations 

(A.7) W' = A p (x,X)W 

indexed by a parameter p, and satisfying exponential convergence condition (|A.2|) uniformly 
in p. Suppose further that 

(A.8) \(A P - A p ± ) - {A - A°±)\ <C\p\e- 9 ^, 6>0 

and 

(A.9) \(A p -A°) ± )\<C\p\. 

Lemma A. 4 ([PZ ]IBHZ] 1. Assuming (|A.2|) and (|A.8|) - (|A.9|) . for \p\ sufficiently small, there 
exist invertible linear transformations P+(x, A) = I + Q p + {x, A) and P_(x, A) = I + @ p _(x, A) 
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defined on x > and x < 0, respectively, analytic in A as functions into L°°[0, ±00), such 
that 

(A.10) \(P P - P°)±| < d\p\e-^ for x ^ 0, 

/or any < < 0, some Ci = Ci(6,9) > 0, and i/ie change of coordinates W =: P±Z 
reduces (|A,7p to iae constant- coefficient limiting systems 

(A.ll) Z' = A^(A)Z for x ^ 0. 

Proof. Applying the conjugating transformation — > (P^) _1 W for the p = equations, 
we may reduce to the case that A is constant, and P? = /, noting that the estimate (|A.8j) 
persists under well-conditioned coordinate changes W = QZ, Q(±oo) = /, transforming to 

'"' <\q^a p -A p ± )-(A°-A ± ))Q- 1 \ + \Q- 1 (A p -A ) ± Q-(A p -A°) ± \, 

where 

(A.13) IQ- 1 ^ - A°) ± Q - (A p - A°)±| = 0(|Q - /|)|(^ - A°)±| = 0( e - e ^)\p\. 
In this case, (|A.8|) becomes just 

|^-Ai| ^Cilple"^, 

and we obtain directly from the conjugation lemma, Lemma lA.H the estimate 

\pP - pQ| = \pP -i\< CCxW~ M 

for x > 0, and similarly for x < 0, verifying the result. □ 

Remark A. 5. In the case A^ = constant, or, equivalently, for which (|A.8[) is replaced by 
\AP-A°\ < Ci\p\e~ e \ x , we find that the change of coordinates W = P P ± Z , P p ± := (P°)± l P^, 
converts (|AT7|) to Z' = A°Z, where P± = I + ® p ± with 

(A.14) [0§.| < CCi\p\e~~ e \ x \ for x ^ 0. 

Taa£ is, we may conjugate not only to constant- coefficient equations, but also to exponen- 
tially convergent variable- coefficient equations, with sharp rate (IA.14D . 

In the general case A p ± ^ A°±, we may still conjugate (|A.7|) to Z' = A°Z 6y the change 
of coordinates W = P±Z , P± := (P°)± 1 Q±P±, where Q± defined by 

(A.15) Q' = A P ± Q - QA° ± 

conjugates the constant- coefficient equation Y' = A p ±Y to X' = A±X, obtaining bounds 

(A.16) Pl = I + 9 p ± , |ejt|<CCi|p| for|x|<C 

valid for finite values of x. However, in general, Q± grow without bound as x — ► ±00. 
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Remark A. 6. As observed in provided that the stable/unstable subspaces of A\/A p _ 

converge to those of A^_/A°_, as typically holds given (|A.9p - in particular, this holds by 
standard matrix perturbation theory JWjj if the stable and unstable eigenvalues of A± are 
spectrally separated- (|A. 10|) gives immediately uniform convergence of the Evans functions 
D p to D° on compact sets of£l, by definition (|A.3|) , 

A. 3 The tracking lemma 

Consider an approximately block-diagonal system 

(A.17) W = j£j (x,p) + 5(x,p)e(x,p)W, 

where is a uniformly bounded matrix, 5(x) scalar, and p a vector of parameters, satisfying 
a pointwise spectral gap condition 

(A. 18) min a(dtMf) - max o-(KAff) > n(x) for all x. 

(Here as usual $IN := (1/2) (./V + N*) denotes the "real", or symmetric part of N.) Then, 
we have the following tracking/reduction lemma of |MaZ3l IPZ] . 

Lemma A. 7 (|MaZ3l IPZ] ) . Consider a system (|A.17|> under the gap assumption (1A.18|) . 
with e uniformly bounded and r\ G L\ oc . If swp(8/r])(x) is sufficiently small, then there 
exist (unique) linear transformations &i(x,p) and <&2( x ,p), possessing the same regularity 
with respect to p as do coefficients Mj and 50, for which the graphs {{Z\, &2Z1)} and 
{(&1Z2, Z2)} are invariant under the flow of (|A.17j) . and satisfy 

(A. 19) sup sup |<J> 2 | < Cswp(8/r]) 

and 

f+CO fX 

(A.20) \§\{x)\<C \ e f y v{z)dz 5(y)dy, |*f (x)\ < C / e** ~ v(z)dz 6(y)dy. 

Jx J— 00 

Proof. By the change of coordinates x ^ x, 5 ^ 5 := 5 / n with dx/dx = r)(x), we may reduce 
to the case 77 = constant = 1 treated in [MaZ3] . Dropping tildes and setting $2 := V^f 1 > 
where (?/4> satisfies ()A.17p . we find after a brief calculation that <3?2 satisfies 

(A.21) $ 2 = (M 2 $ 2 -$ 2 Mi) + <5Q($ 2 ), 

where Q is the quadratic matrix polynomial Q(&) := ©21 + ©22^ — $@n + Q@i2&- Viewed 
as a vector equation, this has the form 

(A.22) $ 2 = M$2 + 5Q($2), 
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with linear operator _M<3? := M2& — $M\. Note that a basis of solutions of the decoupled 
equation <£' = may be obtained as the tensor product = (fxjf of bases of solutions of 
(f)' = M2<j) and <f>' = —M*(f>, whence we obtain from (|A.18|) 



(A.23) e mz < Ce-^, for z > 0, 

or uniform exponentially decay in the forward direction. 

Thus, assuming only that $2 is bounded at —00, we obtain by Duhamel's principle the 
integral fixed-point equation 

(A.24) <D 2 (x) = T<D 2 (x) := f e M ^S(y)Q(^ 2 )(y) dy. 

J —00 

Using (|A.23|) . we find that T is a contraction of order 0{5/rj), hence (IA.24D determines a 
unique solution for 5/rj sufficiently small, which, moreover, is order 5/rj as claimed. Finally, 
substituting Q($) = 0(1 + j$| 2 ) = O(l) in (|AT24j) . we obtain 



|*2(a;)| <C f e^ x -^S(y)dy 
J —00 



in x coordinates, or, in the original ^-coordinates, (|A.20|) . A symmetric argument establishes 
existence of <£>i with the asserted bounds. Regularity with respect to parameters is inherited 
as usual through the fixed-point construction via the Implicit Function Theorem. □ 

Remark A. 8. For r) constant and 5 decaying at exponential rate strictly slower that e^ x 
as x — > +oo ; we find from (|A.20j) that $2(3) decays like 5/rj as x — » +00, while if 5(x) 
merely decays monotonically as x — > — oo ; we find that ^(x) decays like {5/rj) as x — > —00, 
and symmetrically for 

Remark A. 9. 1. A closer look at the proof of Lemma A.l shows that, in the approximately 
block lower-triangular case, 5&21 not necessarily small, there exists a block-triangularizing 
transformation $1 = 0(sup\5/r/\) « 1, under the much less restrictive conditions 

sup (VA7|(I©hI + 16221)) < 1 and sup(|5/r/||e 21 |) « ' 



sup \5/r]\ 



2. Similarly, in the standard, approximately block- diagonal case, an examination of the 
proof shows that bounds (IA.19I) may be sharpened to 
(A.25) 

sup I $1 1 < Csup((5/r/)( sup|6i 2 | + sup(<5/r?) sup(|6n [ + [©22D + sup(<5/?7) 2 sup |0 2 i| 



sup |^> 2 1 < Csup((5/r/)(^sup |9 2 i| + sup(5/rj) sup(|0n| + [©22 1) +sup(5/r]) sup|6i 2 

Remark A. 10. An important observation of fMaZS\, \PZ$ is that hypothesis (|A.18P of 
Lemma \A.7 may be weakened to 

(A.26) mincr^Mf) - max <r(9£Mf ) > t](x) + a{x,p) 
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with no change in the conclusions, for any a satisfying a uniform L 1 bound \a(-,p)\^i < C\. 
(Substitute e Mx < Ce Cl e~ r,z for (|A.23|) . with no other change in the proof.) More generally, 
(|A.27p may be replaced in Lemma A.l\ by 



(A.27) W=P 1 + ai M2 ^ a2 )(x,p) + (5e + /3)(x,p)W, 

where \ctj\j\ are bounded and \f3\ii is sufficiently small, with the conclusion that 
(A.28) sup|$i|, sup|$ 2 | < C(sup(<Vr?) + |/3| L i). 

(T/te additional term f*^ e M ^-^ Q ^ 2 ){y) d y, ■= 021 + /3 22 $ - */?u + $/?i 2 $, 

now; appearing in the righthand side of ()A.24p is contractive for small.) These allow 
us to neglect commutator terms in some of the more delicate applications of tracking: for 
example, the high-frequency analysis of \MaZS$ . or the analysis of case lib in Section \6.2.2[ 
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